The Fokker-Planck operator at a continuous phase transition 
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I consider a physical system described by a continuous field theory and 
enclosed in a large but finite cubical box with periodic boundary conditions. 
The system is assumed to undergo a continuous phase transition at some 
critical point. The q> A theory that is a continuous version of the Ising model is 
such a system but there are many other examples corresponding to higher spin, 
higher symmetry etc. The eigenfunctions of the corresponding Fokker-Planck 
operator can be chosen, of course, to be eigenfunctions of the momentum 
operator. It is shown that the eigenvalues of the FP operator, corresponding to 
each eigenvalue q of the momentum operator, evaluated at a transition point of 
the finite system, accumulate at zero, when the size of the system tends to 
infinity. There are many reasonable ways of defining a critical temperature of 
a finite system, that tends to the critical temperature of the infinite system as 
the size of the system tends to infinity. The accumulation of eigenvalues is 
neither affected by the specific choice of critical temperature of the finite 
system nor by whether the system is below or above its upper critical 
dimension. 



The property of critical slowing down [1,2] is known for a long time from 
experimental [3-7] and numerical work [8-11]. The quantitative description is in 
terms of characteristic decay times or alternatively "characteristic frequencies" co q , 

that govern the decay of a disturbance of wave vector q. It was found that at the 
transition, co q behaves as some positive power of q , co q <x q z . This implies that the 

larger the scale of the disturbance the longer it takes to decay and a divergent scale 
results in a divergent decay time. The evaluation of the exponent z was the subject of 
theoretical work using a number of different approaches but all based on stochastic 
field equations of the Langevin type to describe the dynamics of the system [12-18]. 
More recently, it was suggested that not only do the characteristic frequencies tend to 
zero with q but the decay function itself becomes slower than an exponential (e.g. 
stretched exponential) for long times, co q t » 1 [19-21]. That property was shown to 

hold for quite general Langevin field equations including in addition to critical 
dynamics, equations of the KPZ type. The derivation of stretched exponential decay 
or any other form of slow decay has to rely on approximations. As will become 
evident later (eq.(12)) the time dependent structure factor needs evaluations of the 
eigenstates and eigenvalues of a "quantum field Hamiltonian" (eq.(10)) obtained by a 
standard similarity transformation from the corresponding Fokker-Planck operator. 
The accuracy with which these have to be evaluated has to increase with the time 
argument of the structure factor. Therefore the problem of very long time decay is 
extremely difficult and it is important to have as many exact results as possible in 
support of approximate derivations. From eq. (12) it will become clear also that a 
necessary condition for slow decay of a disturbance of wave vector q is to have 
"enough" eigenvalues corresponding to states carrying momentum q in the vicinity of 
zero. It has been shown in a recent publication [20], that there must exist an 
eigenvector of the FP operator, at the transition point, that carries momentum, q and 
has an eigenvalue that is as close as we wish to zero. In the present article, I go 
beyond that and show that for each q the eigenvalues of the FP operator accumulate 
at zero for the infinite system. The proof I present relies on the simplest established 
lowest order results of finite size scaling equilibrium theory and does not involve 
even the next order corrections such as the dependence of a properly chosen critical 
temperature of the finite system on its size [22, 23]. To my knowledge, those results 



have never been proven rigorously. Nevertheless their use as conjectures in the 
present proof yields information about a non equilibrium problem that is 
considerably more difficult and that information is exact provided that those widely 
accepted equilibrium results are indeed correct. 

Consider a system described in terms of a scalar field q> and enclosed in a cube with 
periodic boundary conditions. The static statistical properties of the system, are 
assumed to be given by the Gibbs distribution,^ cc exp[-W], where W is the 

classical dimensionless Hamiltonian of the system. I have in mind Hamiltonians that 
are even functionals of the field and are generalizations of the q> 4 theory that in term 
of Fourier components of the field has the form 



andQ is the volume of the system. The momentum indices are assumed to be bound 
from above by some high momentum cut off and their allowed values are given by 



where L is the linear size of the system and n a vector of integers. 
The following discussion and the proof I present is not limited, however, to the above 
typical model. The discussion clearly covers other systems with different dependences 
of W on the field (e.g. systems describing higher spins) or systems in which the field 
is not a scalar (e.g. systems with O(n) symmetry). As will become evident , the only 
property that I use is that there is a critical point in the large volume limit and that at 
that point and in that limit the structure factor < <p q <p_ q > or its analog in the case of 

vector fields, diverges at small g's as q raised to a power smaller than -1. 
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The dynamics is assumed to be described in terms of a set of Langevin field equations 
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where the noise rj q obeys 
n (f )) = and (n (f >7_ q ( = 2yS( - ^ ■ (5) 



A standard procedure leads to the Fokker-Planck equation for the probability 
distribution of the fields, P, 
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P = OP, (6) 



where O is the FP linear operator. The distribution at equilibrium, P eq is the eigenstate 
of the Fokker-Planck operator, O, with eigenvalue zero. It follows directly from 
equation (6), that 

P eq oc Q xp[-W] , (7) 
is the Gibbs distribution. 

I will be interested in the time dependent structure factor (q>_ q (0)^ q ( yj , where the 
meaning of the average is as follows: cp_ q is measured at time t=0 at equilibrium, the 
system is then allowed to evolve freely and ^ q is measured at time t. The 
mathematical form of the above statement is 

<P- q (Ok, (0) = iD^DcpY^ V }p?p{<pV ; p<» , I (8) 
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where P eq is normalized and P{<p (2) ;<p (l) , } is the solution of the Fokker-Planck 
equation (5) for a distribution in the variables (p {1) with initial condition 



iV 2) ;^o}=n^ (2) -^ a) )- (9) 
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A standard transformation, P = P^ 2 y/ , which induces a similarity transformation on 

O, brings the Fokker-Planck eq. (6) and the definition (8) into forms more familiar 
from quantum mechanics. 
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y/ = - Hy/, 



(10) 



and 



<P-A0)m( ) = (0 m exp[-!}{t]m 



(11) 



The "Hamiltonian", H, obtained by the similarity transformation from O, is Hermitian 
and all its eigenvalues are positive apart from the eigenvalue zero, which corresponds 
to the ground state |0) . Since the system is translationally invariant, ^commutes with 

the momentum operator so that the eigenstates of 'H can be chosen to carry definite 
momenta. 



Let { | « q ) } be the set of eigenstates of 9f carrying momentum q and [An q } the 
corresponding eigenvalues. It is easily verified that 
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(12) 



(This follows from the fact that m q \0j is an eigenstate of the momentum operator with 
momentum q). The spectrum of ^limited to the states { |w q ^ }, [jin q }, has thus direct 
bearing on the decay of the time dependent structure factor. For example if the 



spectrum \Xn q \ has a gap A q , the decay of the time dependent structure factor is 



faster than e q 
Define next 



^W=Z|(" q KI°)| 2 ^K -^->k(°))- ( 13 ) 



The time dependent correlation function can be written as 

t 

K(oK(0) = W ex pM ], (14) 





It is clear, that for a finite system, it will always be possible to find a region X > 
including the origin and depending on the size of the system, such that p q (X) 

vanishes identically in that region, implying a decay in time that is faster than some 
exponential. This may be characterized by a lifetime that depends on the size of the 
system. The question is how does p q (X) behave in the limit of an infinite system. If 

in that limit p q (X) does not vanish too fast as X tends to zero, the decay will be 

slower than exponential (e.g. stretched exponential, power law etc.) The definition of 
p q (X) implies that in order for it not to vanish too fast as X tends to zero two 

conditions have to be met. First, the eigenvalues Xn q have to accumulate at zero as the 

size of the system tends to infinity. The second condition is that the matrix elements 
in eq.(13) do not vanish too fast as the corresponding eigenvalues tend to zero. The 
actual long time dependence of the decay is strongly affected by the behavior of those 
matrix elements. The accumulation of eigenvalues to be discussed next is thus a 
necessary condition for slow decay but it does not determine the form of the decay. In 
fact, it is not even sufficient. The structure factor in the linear system to be discussed 
in the following as a prelude to the general case decays exponentially in spite of the 
relevant eigenvalues accumulating at zero. The reason is that 
| n (1) q >= q> q | > / < ^ q <p_ q > 1/2 is an exact normalized eigenstate of the Hamiltonian 

so that apart from the matrix element < n (1) q \<p q | > all other matrix elements in eq. 
(12) vanish. The proof I'll present concerns only the first necessary condition that as 



the size of the system tends to infinity the eigenvalues of O accumulate at zero for 
each q. This has no bearing on the actual form of long time decay. 

The strategy of the proof is as follows: (a) The set of many "excitation" states 
carrying momentum q is defined and the corresponding expectation values of the 

"Hamiltonian" .Tf are obtained and shown to accumulate at zero for anyq. (b) The 
relation between those expectation values and the expectation values in an orthogonal 
set of states corresponding to the original "almost orthogonal" set of many 
"excitation" states is outlined, (c) It is shown how the above implies that the 
eigenvalues of the "Hamiltonian" carrying momentum q accumulate at zero. 

Before proceeding to the proof in the general case, I'll describe first a number of 
common results from equilibrium theory that will be used as conjectures in the proof 
and define some useful definitions. Then I'll show how the accumulation of 
eigenvalues of 9f works for the soluble, linear case. This will give some hints and 
motivations for the general non-linear case. 

(1) I keep the value of q fixed (eq. 3) as the size of the system is increased by 
considering only multiplication of the original size, L, by an integer. (2)There are 
many reasonable ways to define the critical temperature of the finite system. It may be 
taken to be the temperature at which the equilibrium structure factor at the smallest 
possible nonzero q is maximal. It may be chosen as the temperature where the specific 
heat or one of its derivatives with respect to temperature is maximal. For our purpose 
it is enough even to choose it as the transition temperature of the infinite system. (4) 
Regardless of the choice of the finite size critical point the structure factor has the 
scaling form. 

(<p q <p_ q ) = Aq-^ fyrqL»l, (15) 

to leading order in L and for q small compared to the cut-off momentum q . Note that 
since q is kept fixed the condition qL »1 can always be attained. Note that that is 
true also above the upper critical dimension where 77=0 regardless of the fact that 
hyper scaling does not hold. Higher order correlation functions will also be needed. I 
will need (cp^ ...<p k ^ , where m is even and = but no subset of the k 's sums 



up to zero. For a large finite system, these correlations have at the critical point the 
scaling form 



(<Pk-<Pk)\ = A (™) l4 L )/2 /(k;,-,kj, (16) 

to leading order in L. For |k|'s small compared toq Q , that obey also 
|k|Z »1 the function /is homogeneous. Namely, if we scale all its variables by a, the 
function of the scaled variables is related to the function of the original variables by 
/(«k 1 ,...,«k m ) = «- m ( 1 -" /2 V(k 1 ,..,kJ. (17) 

The pre factor A(m) is of combinatorial origin and we will need to assume that there 
exists a finite b such that \nA(m) < (bm)hi(bm) . 

(4) A set of momenta, S k = that can contain, in principle, the same 

momenta a number of times, is said to be irreducible if it does not contain real 
subsets, the momenta of which sum up to zero. 

(5) A set S k conjugate to 5^ is the set {-1 1V ..,-1 A }. 

(6) The sum S l © S 2 of the sets S t and S 2 , contains all momenta appearing in 
5, and S 2 . The number of times each momentum appears in the sum is the sum of the 

number of times it appears in each set separately. 

(7) The many "excitation" set is the set of states of the form 



m 

n( l i—iJ=n^l°)- vi/ q i o >' (is) 

m 

where ^1. =q and where the set {lj,...,l m } 



i=l 



is irreducible. 

A "physical proof would probably proceed along the following lines. First the 1,. "s 
are chosen to equal q/m . Then eqs.(21) and (22) are used to obtain the corresponding 
expectation value, q^ 2 ^ l{Am (l ~ ll) ) . The third step would be to assert that because of 
eq.(17) ,the "many excitation" states form an orthogonal set in the large volume limit 



(that is also orthogonal to the ground state at any finite size , because it carries a 
different momentum), the eigenvalues of the "Hamiltonian" {2« q } accumulate at zero. 

The "mathematical proof proceeds basically along the same lines but has to take into 
account two difficulties. The first is that if q is an allowed momentum, q/m is 
generically not allowed. The second is much more serious. The assertion that the set 
of representative states is orthogonal in the large volume limit is problematic. Many 
examples exist in which the fact that the scalar products of normalized states tend to 
zero with the size of the system does not imply consequences expected from 
orthogonal sets in finite systems. 

The proof for the Gaussian case (w=0) is straight forward. In that case, the states 
(18) are exact eigenstates of the "Hamiltonian" .Tfwith eigenvalues. 

r q (i l ,...,\ m )=(r/2)Y J if, (19) 

In the following I will introduce certain restrictions on the size of m , relating it to q .1 
will use these restrictions for the Gaussian case as well for the general one. These do 
not limit the final result about the accumulation of the eigenvalues but enable an exact 
proof under the conjectures about equilibrium correlation functions at the transition. 

For any given q the size of the system can be chosen large enough so that n, the 
absolute value of the vector n corresponding to q, is large enough so that after 
choosing a number < 6 < min(l/2,# ) , where O will be specified later (eq. (25)), 

the integers m in the range (n e ,\.5n e ) are large enough , so that the following will 
hold: (a) It is possible to choose for each m a representative state of the form (18) such 
that the sizes of the vectors l ( . m defining that representative are bound from above by 
^/(w-l).(The idea is to choose all the l f 's as close as possible to q/m since the 
choice 1. = q/m for all i yields the lowest possible value for the eigenvalue in the 

sector of "m excitations" states. The above obvious choice is impossible, however, 
because if q is an allowed momentum q/m is generically not allowed. The 

conditions on m ensure that it is possible to choose thel,- 's in such a way that the 
corresponding eigenvalue (eq.(19)) departs in a controlled way from the absolute 
minimum of the expression ).b) Furthermore, the set Sm x © Sm 2 is irreducible for any 



m j and m 2 in the range. (Note that this is not required for the proof in the Gaussian 
case but will simplify the proof in the general case.) Consequently, the representative 
eigenvalues are bound from above 



f (l 1 ,...,\ m )<(r/2)-^q 2 <( r /2)^, to leading order in n. (20) 

\m-\) n 



Since n = qLI2n , the meaning of eq. (20) is that there are at least ^(qL/27r) d states 

2 

with eigenvalues smaller than {y 1 2) - — - — — . Namely, as the size of the system 

(qL 1 2n) 

tends to infinity, the eigenvalues of H accumulate at zero for any q . 

The proof for the general case will proceed in two steps. The first is to show that the 
expectation values of H in the many "excitations" states (18) accumulate at zero. The 
second is to show that this accumulation implies that the eigenvalues of H also 
accumulates, at zero for any q when the size of the system tends to infinity. 



The representative states are exactly those chosen before for the Gaussian case. The 
corresponding expectation values of the Have 
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The above is an exact expression that holds for a system of any size. The main points 
that lead to the first equality above are invariance under reflection 
«0 | W_ q JC ^ q | 0> = <0 | W q H W_ q | 0) and H |0) =0. The second equality in (21) 

d 2 

above is due to the fact that the double commutator [ v P_„,[y' ^q]] equals 

d(p x d(p_ x 

„ d*¥„ d*¥ . 

2> .It may seem that the last expression on the left hand side of (21) does 



not depend on the non-linear coupling u in the "Hamiltonian". This is not true of 
course, because the ground state |0) depends on u. To leading order in the size of the 
system 

m m 

(° I n <p- h <p h I °> = (° k <p- h I °>(° I n n Pi, I °> > ( 22 ) 

i'=l i* j 



The corrections are terms of the same order of magnitude multiplied by a 
combinatorial factor of the order of m and divided by the volume of the system, the 
correction is of the relative order of L w ~ d , that tends to zero for any d (the 
dimensionality of the system) as the size of the system tends to infinity. Thus, for a 
large enough system, the approximation (22) can be used to yield 



\n J 



(23) 



(Note that qL here is still large in spite of the large m. In fact it is 2m ( ' and 
therefore much larger than 1 so that eq.(15) is used for the structure factor.) This 

means, that there are -^(qL/2^) states, at least, for which the expectation values of 

H&xq smaller than— (3I2) 9 — j—^- = s(q,L) . Since ri is always smaller than 

1 this implies that as the size of the system tends to infinity, the expectation values of 
^"accumulate at zero for any given q. Note that so far I have used only the scaling 
form of the two point function (eq.(15)) as the only input conjecture in the proof. 

If the representative states above were orthogonal, it would be easy to complete the 
proof that the eigenvalues of H corresponding to states carrying momentum q 
accumulate at zero. This, however, is not the case. The purpose of the next part of 
the proof is to show that the set of representative states, described above, can be used 
to construct an orthogonal set of states, by a Grahm-Schmidt orthoganolization 
procedure, such that the expectation values of H in the states of that set are identical to 
leading order in the size of the system to the corresponding expectation values 



calculated above. Because of the irreducibility of any set, Sm 1 ® Smi the scalar 
products of pairs of the representative states are given by eq. (16). Denote, for 
simplicity the many "excitation" representative states by y/ m and the smallest and 

largest w's between n e and \.5n 9 by m s and m l respectively. Further denote by 

max and min the maximal and minimal values respectively of |(^ ; , y/. )| for i ^ j in the 

range m s < i j < m, . 

The orthogonal set {% m } is obtained by constructing linear combinations of the y/'s 

i-i 

Z i =Vi + T j ayVj , (24) 

j=m s 



It is not difficult to show that if min can be made as large as we wish compared to 
[max] 2 n e I A"" i! 2 ^"" when the size of the system is increased, a .,. is given to leading 

order in the size of the system by --, -r for j>i. Use of equations (15) and (16) 

enables to show that the inequality relating min and max can be made to be obeyed 
for large enough systems. It is easy to show now that the required inequality relating 
max and min is obeyed if 

0<(d-2 + r])/2b = o . (25) 

This choice of the range of 6 ensures that for j in the range, [(^^J is 

monotonically decreasing as a function ofz' + j , for large enough L. It also ensures 
that, 

(Zi ,Zi)= (Wi , ¥i ) and (% t \ 9f Zt)= (Wi ,Hy/), (26) 



where relative corrections are at most of the order of [max] 2 7/ / min 2 A" L (2 n)n 



(In fact, it may seem that a weaker condition on 6 , 9 < (d - 1 + r/ 1 2)1 2b ( r/ is always 

smaller than 2), suffices, to obtain eq. (26) but since the monotonicity of |(^,^ y )| as a 

function of (/' + j) simplifies matters considerably, the somewhat more stringent 
limitation is used.) This implies that to leading order in L, the expectation values of H 
in the orthogonal set {% m } are the same as the expectation values in the many 

"excitation" states. Note that in the proof of the existence of the orthogonal set that is 
"close" to the non orthogonal set of many excitations I used the input scaling form of 
the many point correlation functions (eq.(17)) . The final stage of the proof is to 
show that the accumulation of the eigenvalues of the "Hamiltonian" follows. Assume 
that this is not the case. If so since the eigenvalues of 9fare positive, there exists a 
finite gap ,A q , independent on the size of the system (for large enough L) and there 

are only a finite number M (independent on the size of the system) of eigenstates of 
H with energies in the gap that tend to zero with the size of the system . Denote these 
normalized eigenstates by O.and the corresponding eigenvalue bX q and any 

normalized vector vector belonging to the subspace spanned by eigenstates of the 
"Hamiltonian" with energies larger than the gap by O a . Consider next the states 
{Xi } with expectation values of the Hamiltonian smaller than s{q,L) .1 will express 
each of these states by 

X=Ifyl>,+W- (27) 
The expectation value corresponding to Xi is bounded by 



On the other hand J u i < s(q,L) . The conclusion is that| b i | 2 < £ ^'^ , Therefore, the 

projections of the states % i on the M dimensional space spanned by the O, 's are 
almost normalized and almost orthogonal to each other. The last result is impossible 



of course in view of the enormous number of sates % t with expectation values less 
than s (q, L). This completes the proof. 

The proof given here concerns only the accumulation of eigenvalues of the FP 
operator at zero in the thermodynamic limit. This result is exact and general under the 
conjecture that the established results for equilibrium correlation functions at the 
transition are correct. It is thus a proof of the existence of a necessary condition for 
slower than exponential decay at the transition point. It has no bearing however on the 
actual form of decay whether exponential, stretched exponential or power law. The 
form of the decay depends on certain matrix elements that are problem specific. I 
hope to discuss this in future work. An interesting question that has been raised a 
number of times in connection with the above description is whether it is relevant to 
the glass transition. The first tendency is to say that it is irrelevant, because in the 
glassy state translational symmetry is broken and that was used heavily in the above. 
Furthermore, the glassy phase is supposed not to be in equilibrium and the above 
treatment is based on equilibrium dynamics. It seems, however, that there is some 
chance of success by considering the liquid phase close to the transition. I hope to 
come back to this in the near future. 
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